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Abstract 

We study the potential impact of detecting the inflationary gravitational wave background by 
the future space-based gravitational wave detectors, such as DECIGO and BBO. The signal-to- 
noise ratio of each experiment is calculated for chaotic/natural/hybrid inflation models by using 
the precise predictions of the gravitational wave spectrum based on numerical calculations. We in- 
vestigate the dependence of each inflation model on the reheating temperature which influences the 
amplitude and shape of the spectrum, and find that the gravitational waves could be detected for 
chaotic/natural inflation models with high reheating temperature. From the detection of the grav- 
itational waves, a lower bound on the reheating temperature could be obtained. The implications 
of this lower bound on the reheating temperature for particle physics are also discussed. 
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I. INTRODUCTION 



Gravitational wave, which is the gravitational counterpart of electromagnetic wave, has 
eluded detection. Since gravitational waves interact very weakly with matter, the Universe 
viewed with the gravitational waves appears much transparent than that with photons. 
Therefore, the detection of the gravitational waves gives us a snapshot of the very early 
Universe. Since we now already know that the cosmic microwave background (CMB) spec- 
trum has been a powerful tool for cosmology and the analysis of it has revealed the various 
components of the Universe, it is not difficult to imagine that another new era of cosmology 
will come after the discovery of gravitational waves. At the time of the discovery of CMB 
in 1965, who could have imagined the present situation of cosmology. It is important to 
prepare for the coming age of another precision cosmology. 

Among the processes that occurred in the very early Universe, inflation (including re- 
heating) is the most important epoch. Inflation was proposed as the most natural solution 
to the difficulties of the standard big-bang cosmology, such as the horizon problem and the 
flatness problem It also generates primordial gravitational waves (tensor modes) ^ 
as well as the primordial density fluctuations (scalar modes) {3]. The latter has already 
been observed as the cosmic microwave background anisotropics by the Cosmic Background 
Explorer satellite {4] and by the Wilkinson Microwave Anisotropy Probe (WMAP) satellite 
5! . The detection of the gravitational waves generated during inflation could determine the 
energy scale of inflation directly, which is essential to explore the physics behind inflation 

In Ref. [7 1, we calculated for the first time the most precise power spectrum of the 
inflationary gravitational wave background for several inflation models with perturbative 
reheating. Our numerical approach has enabled us to obtain the precise amplitude of the 
spectrum over all frequencies. We compute the spectrum all through the evolution of the 
Hubble expansion rate, which determines the amplitude of the gravitational waves, by fol- 
lowing the dynamics of the inflaton scalar field and its decay into radiation (reheating). We 
fully take into account of the changes of the effective number of degrees of freedom during 
;he radiation-dominated era which cause the suppression of the spectrum at high frequencies 
8j. Therefore, all factors which affect the amplitude of the gravitational wave background 
spectrum are precisely reflected in our calculation. 
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In this paper, we aim to forecast the future prospects of the direct detection of the 
gravitational wave background by next generation satellite detectors, such as the DECi-hertz 
Interferometer Gravitational wave Observatory (DECIGO) Q, Q and Big-Bang Observer 
(BBO) 11]. While the measurement of CMB B-mode polarization, which is the indirect 
signature of primordial gravitational waves, probes gravitational waves at the present horizon 
scale, the direct detection observes gravitational waves at a much smaller scale, which would 
contain unique information about the very early Universe. In particular, we should stress 
that the shape of the gravitational wave spectrum around the target frequency (~ 0.1Hz) is 



sensitive to the physics of reheating 



121 ] . In Ref. 13j . the detectability of the gravitational 



wave background is estimated taking into account the dependence of the spectrum shape 
on reheating temperature. Here, we reexamine the detectability by calculating the signal- 
to-noise ratio (SNR) with the precise prediction of the spectrum amplitude and the specific 
noise spectra of DECIGO and BBO. We study the dependence of the spectrum on reheating 
temperature as well as on models of inflation. In addition, we discuss the implications of a 
possible lower bound of the reheating temperature obtained from the future direct detection 
by DECIGO/BBO for particle physics. 

Note that our investigation is carried out assuming reheating via perturbative decay of 



the inflaton fie 
preheating 14 



d, namely we do not consider nonperturbative effects during reheating, called 



151 ] . If the nonperturbative effects are dominant in the reheating process, 



the picture of reheating and its effect on the gravitational wave background are significantly 



different from those of perturbative reheating 1 6M 1 31] . In this paper, since we would like to 



provide a conservative estimate of the gravitational wave background and its detectability, 
we only consider perturbative processes which always exist and are directly related with the 
reheating temperature. 

The outline of this paper is as follows. In Sec. [Ill the signal-to-noise ratio expected in 
future experiments is estimated for four inflation models; chaotic inflation with quadratic 
and quartic potentials, natural inflation and hybrid inflation. First of all, we present the 
method to calculate the signal-to-noise ratio and the current specification design of DECIGO 
and BBO. Then the procedure of the numerical calculation is described briefly, which is used 
to obtain the amplitude of the spectrum. In Sec. 1111} we calculate the signal-to-noise ratio 
using the detailed experimental specification and the precise prediction of the spectrum 
amplitude. The implications of the lower limit of the reheating temperature on particle 
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physics are also discussed in relation with the gravitino problem in Sec. IIVI A summary is 
given in Sec. [V] 



II. ESTIMATION METHOD OF THE SIGNAL-TO-NOISE RATIO 

A. Correlation analysis for detection of a stochastic gravitational wave background 

Cosmological gravitational waves are described as tensor perturbations in the Friedmann- 
Robert son- Walker metric as ds 2 = —dt 2 + a 2 {t)(5ij + hij)dx % dx^ , where a(t) is the scale factor 
of the Universe. The tensor perturbation hy can be expanded into its Fourier components 

as 

hAt, x) = £ J ^4( k )^(*K k ' x , (i) 

where the polarization tensors ep x satisfy symmetric and transverse-traceless conditions 
and are normalized as £\ . e^(e^)* = 25 xx ' . The intensity of a stochastic gravitational wave 
background is characterized by the dimensionless quantity, GW = (dp GW /d\iak)/ p c , where 
the critical density of the Universe is defined as p c = 3H 2 /8ttG with the Hubble expansion 
rate, H = (da/dt)/a. The energy density of the gravitational waves pew is given from the 
00-component of the stress-energy tensor as pew = ((dthij) 2 + {Vhij / a) 2 ) / {QAnG) . Then 



""4 s sEi* < 2 > 



Hqw can be expressed in terms of the Fourier component as 19] 

2 If 

12 \aH J 7T 

In future satellite missions like DECIGO and BBO, the analysis of a stochastic gravita- 
tional wave background would be performed by taking the cross correlation of the outputs of 



gravitational wave detectors |2CH22j]. The signal-to-noise ratio for the correlation analysis is 
given in terms of an expected (theoretical) form of Oqw(/)) and the functions related to the 
experiment design, such as the noise spectrum Si j(f) and the overlap reduction function 



lutf) as [23| 

w=<|)V.Ef,» ( 3) 

where / = k/2n is the frequency of gravitational waves, H is the present Hubble expansion 
rate, T Q b s is the duration of the observation time. The subscripts I and J refer to independent 
signals obtained at each detector, or observables generated by combining the detector signals. 



DECIGO is planned to be a Fabry-Perot Michelson interferometer with an arm length of 
L — 1.0 x 10 3 km 0, Q. In this case, the SNR is calculated with the noise spectral density 



of the two interferometers, which are assumed to be the same and given by [20, 

Sl(f) — 82(f) — S s hot + S'acccl + S ra _d, (4) 

where the shot noise is given as S s hot = 5.29 x 10~ 42 (1 + / 2 // c 2 )(L/km)~ 2 Hz _1 , the ac- 
celeration noise is S^ccei = 4.0 x 10 _46 (//Hz) _4 (L/km) _2 Hz _1 and the radiation pressure 
noise is 5 rad = 3.6 x 10~ 51 (//Hz)~ 4 (l + / 2 // 2 )- 1 Hz' 1 . 1 The cutoff frequency is given by 
f e = 1/(4FL) with the fineness for the DECIGO detector, 7 = 10. 

On the other hand, BBO would adopt a technique called time-delay interferometry, in 
which new variables (I = A, E, T) are constructed to cancel the laser frequency noise. The 



noise transfer functions for the time-delay interferometry variables are given as 24], |25] 

S A (f) = S E (f) = 8sin 2 (//2)[(2 + cos/Vahot 

+2(3 + 2 cos / + cos(2/))S accel ], (5) 

S T (f) =2[l + 2cos/] 2 [S shot + 4sin 2 (/72)S accel ], (6) 

where / = 2irLf. In the case of BBO, the arm length is L = 5.0 x 10 4 km, 
and the noise functions are S s hot = 2.0 x 10~ 40 /(L/km) 2 Hz _1 and S^cei = 9.0 x 
10- 4 7(27r//Hz) 4 /(2L/km) 2 Hz- 1 . 

The overlap reduction function r )u(f) can be calculated by taking into account of the 
relative locations and orientations of the detectors. We use the results of Ref. [3] for 
FP-DECIGO, and of Ref. U for BBO. 



B. Numerical calculation for the spectrum of the gravitational wave background 

We briefly describe the method we used to obtain the theoretical prediction for the ampli- 
tude of the spectrum Oqw in our previous work (For details, see Ref. 3]). We numerically 
solve the evolution equation for gravitational waves, which is derived from the perturbed 
Einstein equation under the Friedmann-Robertson- Walker metric, 

k 2 

h£ + 3Hh£ + -^ht = 0, (7) 



1 A major improvement of the sensitivity is under consideration. 
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where the over dot describes the time derivative. The initial condition is randomly taken 
from the Bunch-Davis vacuum, of which variance is given as 

where mpi = G^ 1 ^ 2 denotes the Planck mass. 

The important point of our evaluation of the spectrum is that we compute the evolution 
of the gravitational waves with following all through the history of the cosmic expansion 
from inflation to the present epoch. It allows us to calculate the spectrum with no use of the 
slow-roll approximation, which overestimates the amplitude of the spectrum at the direct 
detection scale in some inflation models. The spectrum, which is Taylor-expanded around 
the CMB scale, can overestimate the amplitude by 10 — 20% [3]. The numerical approach 
also enables us to precisely evaluate the effect of the changes in the relativistic degrees of 
freedom during the radiation-dominated era. 

If we assume reheating is proceeded by perturbative decay of the inflaton field into light 
fermions, all the processes from inflation to the end of reheating can be calculated by simul- 
taneously solving the following equations, 

0+(3iJ + r)0 + ^M = O, (9) 
p r + AHp r = T P<t> , (10) 

H2 = i^r^ + Prl (11) 

where V is the decay rate of the scalar field into radiation, V(<f>) is the potential of the scalar 
field, p r is the energy density of the radiation, and the energy density of the scalar field is 
given as p^ = 2 /2 + V ((/)). 

During inflation, the Hubble expansion rate is determined by the dynamics of the scalar 
field which drives inflation. The Universe enters a reheating phase after inflation, and the 
scalar field oscillates at the bottom of the potential decaying into radiation. During this 
phase, the Universe evolves like a matter-dominated Universe in most of the inflation models, 
and turns into a radiation-dominated era after it ends. The exception is, for example, the 
case where the potential has a A0 4 shape at its bottom. In this case, the Universe behaves 
as a radiation-dominated Universe and connects to the subsequent radiation-dominated era 
with no change in the Hubble expansion rate. 
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If the matter-dominated reheating phase exists before the radiation-dominated era, the 



inflationary gravitational wave spectrum is suppressed at high frequencies [12|, [13|. The 
suppression is seen on the modes which enter the horizon during reheating, since the matter- 
dominated phase induces frequency dependence of f~ 2 on the spectrum while a radiation- 
dominated era gives a flat spectrum oc f°. The characteristic frequency, where the change 
of the frequency dependence from f~ 2 to/ arises, is given in terms of the temperature of 



27] 



the Universe at the end of reheating as 

/rh - 0.3 ( ^ RH ) f 9 -^) V2 " 1/3 Hz, (12) 

JR \WGeV J V 220 / V 220 / ' v ; 

where g* and g* s are the effective number of relativistic degrees of freedom contributing to 
the radiation density and the entropy density. The subscript "RH" denotes the value at the 
end of reheating. We take <?*,rh to be ~ 220, which includes degrees of particles in minimal 
supersymmetric standard model. Since both g* and g* s take the same value before electron- 
positron annihilation, ~ O.lMeV, where the temperature of neutrinos becomes lower than 
that of photons, (?* Si rh is also taken to be 220. In the case of^the perturbative decay, the 
reheating temperature T RH is related to the decay rate Y as 



Tffl ^^ 4) (S) 1/4(mpir)1/2 (13) 
Note that the above expressions are valid only for the case where the mass of the in- 
teracting fermion is much smaller than the mass of the inflaton field and also the coupling 
constant is sufficiently small. Otherwise, parametric resonance between the two interacting 
fields may give nonperturbative creation of particles and dramatically shorten the time scale 
of reheating. We do not deal with such nonperturbative particle production during reheat- 
ing, called preheating [l^] . If the nonperturbative growth is sufficient to reheat the Universe, 
the f~ 2 dependence would not arise on the inflationary gravitational wave background and 
it may rather be important to look at gravitational waves origina ting from large inhomo- 



geneities in the matter field during the nonperturbative stage |16l-ll8j. Such gravitational 
waves typically have very high frequency beyond the sensitivity bandwidth of DECIGO and 
BBO, but some models predict an infra-red tail which overlaps or exceeds the inflationary 
gravitational wave background at detectable frequency. However, in this paper, our focus 
is on the effect of reheating on the inflationary gravitational wave background and consider 
only the case of reheating with perturbative fermionic decay. For reheating via fermionic 



decay, the rapid particle production is suppressed by Pauli blocking, while for reheating 
via bosonic decay could be easily completed by the exponential growth of the number of 
particles. Alth oug h the effect of the nonperturbative decay can still be important for the 
fermionic case [15| and could change the time scale of reheating, our investigation, after all, 
turns out to be only the weak coupling case, where the nonperturbative effects are negligibly 
small. The detailed analysis of the effects of nonperturbative processes like preheating on 
the inflationary gravitational wave background will be left to a future work. 

After reheating ends, the Universe enters a radiation-dominated phase. Including the 
effects of g*{T) and g* s (T), the Hubble expansion rate can be written as 29] 

-4/3 



Hi 



9*(T)\ ( 9*s(T) 



g*,o 



g*s,o 



Vt r a 



VL m a + Vtt 



(14) 



where denotes values at the present time, which are given by summing the contributions 
from photons and neutrinos: g* t Q = 3.36 and g* Sj o = 3.90. The change of the Hubble expan- 
sion rate affects the evolution of the inflationary gravitational waves and causes damping at 
the frequency where the direct detection experiments are targeting. The suppression is about 
(g*(T = 10 7 GeV)/^, )(^(T = 10 7 GeV)/^ s , )~ 4/3 = (220/3.36)(220/3.90)" 4 / 3 ~ 0.3. 

We calculate the spectrum of the gravitational waves by numerically solving the evolution 
of each mode with Eq. (J7|). At the same time, we follow the evolution of the Hubble 
expansion, which is calculated by Eqs. ( ITOj) and flTTJ during the inflation and reheating 
phase, and by Eq. ( JT4"j) after reheating. In this paper, the energy density of radiation is 
taken to be Q r h 2 = 4.15 x 10" 5 and the other cosmological parameters are given by the 
maximum likelihood values from the combined constraints of the WMAP 7 yr, BAO, and 
supernova data30[: matter density fl m h 2 = 0.1344, cosmological constant density Q\ = 
0.728, amplitude of curvature perturbations = 2.45 x 10~ 9 , and the Hubble parameter 
h = 0.702. 



III. PREDICTIONS FOR THE DETECTABILITY IN FUTURE EXPERIMENTS 
A. Signal-to-Noise Ratio 

Now we forecast the detectability of the inflationary gravitational wave background with 
DECIGO/BBO for several inflation models: chaotic inflation with a quadratic and quartic 
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potential, natural inflation 3l|, |32| , and hybrid inflation 33] . The potentials are respectively 
given by, 

quadratic : V = -m 2 2 , (15) 



2 

quartic : V = ^0 4 , (16) 



natural : V = A 



1 ± cos 



N(p 

T 



(17) 



hybrid: V = jr(M 2 - Aa 2 ) 2 + -m 2 2 + -g 2 2 cr 2 . (18) 
4A 2 2 

The spectra of these four models and the sensitivity curves of the DECIGO and BBO 
experiments are shown in Fig. [TJ The normalization of the scalar perturbations = 
2.45 x 10~ 9 fixes the value of the potential parameters. It gives m = 1.72 x 10 13 GeV for the 
quadratic potential, A = 1.54 x 10~ 13 GeV for the quartic potential, A = 2.04 x 10 16 GeV for 
the natural inflation model with N — 1 and / = 2mpi, A = 1.33 x 10 16 GeV for the natural 
inflation model with N — 1 and / = mpi, and M = 1.45 x 10 16 GeV for the hybrid inflation 
model with X = 1, g = 8x 10~ 4 , and m = 2.5 x 10~ 7 mpi. 2 The parameter values are searched 
numerically by adjusting the resulting Hubble expansion rate to be H = 100/ikm/s/Mpc. 
Note that these parameters also slightly depend on the reheating temperature, since it is 
related to the length of inflation. The above values are obtained when Trh = 10 7 GeV. 

In Fig. [21 we show the SNR for cross-correlation analysis expected with a 10-year obser- 
vation by DECIGO and BBO. The SNR decreases significantly when the reheating temper- 
ature Trh is lower than 10 7 GeV, since the suppression due to the presence of the reheating 
phase reduces the amplitude of the spectrum at the target frequency of the direct detection 
experiments. The only exception is the case of the quartic potential, in which the Hubble 
expansion rate behaves as a radiation-dominated Universe during the reheating phase and 
the suppression does not arise. 

In Table [H we present the relation between the value of ^gw at the detection frequency, 
/ = 0.2Hz, and the tensor-to-scalar ratio, which is usually evaluated at the CMB scale 
fc CMB = 0.002Mpc as 

r s ^mb) ^ 16e _ 



2 Although there is much freedom in the choice of the parameters for hybrid inflation, here we take one 
example where the parameter M corresponds to grand unification scale. 
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FIG. 1: Spectra of the gravitational wave background for different inflation models, shown with 
the sensitivity curves of DECIGO (dotted) and BBO (solid). The spectra are calculated assuming 
Trh = 10 7 GeV. The cases of Trr = 10 6 GeV and 10 8 GeV are also plotted assuming the quadratic 
potential model. Note that the spectrum lines mean the time-averaged value of Oqw- 



where the slow-roll parameter is defined as e = rripj (16Tc)(y f /V) 2 \k CMB=a H- The reheating 
temperature is set as T RII = 10 9 GeV which is so high that the suppression does not arise 
at the detection frequency. As is clear from the comparison between m 2 2 and A0 4 model, 
the amplitude of the gravitational wave at the direct detection scale is not proportional to 
the tensor-to-scalar ratio r because the tilt of the spectrum ~ — 2e and the higher order 
;erms of the Taylor-expansion becomes important when connecting the two different scales 



34(. This is prominent in models which predict larger r. 

For DECIGO, the inflationary gravitational background could be detected with SNR> 3 
if the inflation model is chaotic inflation with Trh ^ 10 7 GeV. For BBO, the inflationary 
gravitational background could be detected with SNR> 5 if the inflation model is chaotic 
inflation with Trh ^ 2 x 10 6 GeV or natural inflation with / > mp\ and Trh > 10 7 GeV. 
Therefore, from the contraposition, if DECIGO does not detect the inflationary gravita- 
tional wave background, the chaotic inflation model will be excluded unless the reheating 
temperature is lower than 10 7 GeV. The same argument holds for BBO except that it can 
apply to natural inflation. 
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FIG. 2: Signal-to-noise ratio vs reheating temperature calculated for the four different inflation 
models. The upper panel shows the SNR for DECIGO and the lower panel is for BBO. The gray 
region shows the la uncertainty in the normalization = (2.44llQQg2) x 10 -9 from WMAP 7 
yr, which is used to determine the energy scale of inflation. 



Model 




r 






SNR (BBO) 


SNR (DECIGO) 


m 2 4> 2 




0.144 


1.12 x 10" 


-16 


13.3 


3.02 






0.262 


8.75 x 10" 


-IV 


10.5 


2.40 


Natural (/ = 


2mpi) 


0.108 


1.01 x 10" 


-16 


12.0 


2.72 


Natural (/ = 


mpi) 


0.0406 


5.70 x 10" 


-17 


6.79 


1.53 


Hybrid 




0.0104 


2.44 x 10" 


-17 


2.90 


0.646 



TABLE I: The tensor-to-scalar ratio r, the amplitude of the gravitational wave f^cw a ^ / = 
0.2Hz, and the signal-to-noise ratio in DECIGO and BBO for each inflation model. The reheating 
temperature is set as Trh = 10 9 GeV. 
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FIG. 3: Parameter dependence of the signal-to-noise ratio for the natural inflation model. 
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FIG. 4: Parameter dependencies of the signal-to-noise ratio for the hybrid inflation model. The 
stars show the fiducial point, at which parameters are set as A = 1, g = 8 X 10 -4 , and m = 
2.5 x 10 -7 mpi. The dotted line in the top panel shows y = 1. The shaded region means a ^ 
when the observable scale (the present Hubble horizon) exits the horizon. 
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B. Parameter Dependence 



We also investigate the parameter dependence of the SNR, shown in Figs. |3] and [4] by 
setting the reheating temperature to be Trh = 10 9 GeV. For the natural inflation model, the 
amplitude of the gravitational waves becomes larger as / increase. This can be interpreted 
as follows: Since we use the normalization of the scalar perturbations A 2 ^ = 2.45 x 1CT 9 , 
the value of e determines the amplitude of the gravitational waves at the CMB scale, as we 
know A 2 h = 16eA| from Eq. ([19]). 3 



In the case of the natural inflation model with N = 1, e is given as 35j 



1 



167T 





sin(x) 


V / J 


1 + cos(a;) 



(20) 



where we define x = 4>J f. The field value can be written in terms of the e-folding number, 
Af = ln(a end /a), as 



sin 



(I) 



sin ^— — j exp 



m 2 Pl 



(21) 



16tt/ 2 

The end of inflation x en d is defined at the point where the slow-roll condition is violated, 
e = 1, which gives 

1 — 167r(f/mpi) 2 , , 

cos( ^> = (22) 

Combining Eqs. fl20D (J2U and (J22K, we find that e is written as e = l/[(l+x) exp(2jV/x)-x] 
where x = 167r/ 2 /m P1 . This is an increase function of Xi which implies that Oqw increases 
as / increases. 

In the standard picture of hybrid inflation that the field determines the evolution of 
the Universe during inflation, the potential can be recast as V = A[l + (0//i) 2 ] for the 
inflation stage, where we define A = M 2 /(4A) and /x = M 2 /(m\/2X) 37]. Then the slow-roll 
parameter is given as 

1 ' mpi \ v 

1 n ■' (23) 



4tt ^ fi ) [1 + y 2 } 2 
The field value is given in terms of the e-folding number [38| 



\ 



W o{ Vend eX P 



2 N 

Vend + 



mpi 



(24) 



3 Although this relation does not hold at scales well below the CMB scale, the amplitude at the direct 
detection scale is a monotonically increasing function of e for the natural and hybrid inflation models, 
which predict relatively small slow-roll parameters. 
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where we define y = and Wq(x) is the principal branch of the Lambert function, which 
satisfies x = Wo(x)e w °( x \ In the case of hybrid inflation, inflation ends with the waterfall 
field a rolling down when enc j = M/g, which gives 



my2X 

</end = ~ Mr - (25) 

The potential looks more like quadratic with increasing m, which is the mass of the </> 
field. When y > 1, the inflation dynamics becomes the same as the case of the quadratic 
potential. So, we focus our interest on the case of y < 1. The value of m affects both \i and 
y, which respectively decreases and increases with increasing m. From Eqs. ( 123"j) (JMD an d 
fl25]), we find e becomes larger with decreasing /i and increasing y for < y < 1. Hence flew 
increases as m increases as seen in the top panel of Fig. HI 

The parameter A determines the potential minimum of the a field, a m m = M/ Since 
observables are determined by the dynamics of the field in the standard hybrid inflation, 
the value of A has little effect on the amplitude of the gravitational waves as shown in the 
middle panel of Fig. |H 

The parameter g determines the value of en d = M/g, where the waterfall field a starts 
to roll down. The value of y is affected through y e nd = 0end//^ which becomes smaller as g 
increases. Thus, the slow-roll parameter e, which is an increase function of y for < y < 1, 
becomes smaller with increasing g. Hence, SNR oc ^gw decreases as g becomes larger in 
Fig. H 

For too large A or small g, the condition rn 2 /(f)l nd = m 2 M 2 /g 2 M 4 /A {33] is not 
satisfied, which results in that the waterfall field a starts to roll down slowly before inflation 
is dominated by the vacuum energy of the <fi field. The shaded region in the figure means 
a ^ when the observable scale (the present Hubble horizon) exits the horizon, which is 
not the standard behavior of the hybrid inflation model. 



IV. IMPLICATIONS OF THE LOWER LIMIT OF THE REHEATING TEMPER- 
ATURE FOR PARTICLE PHYSICS 

In previous sections, we found that if inflation is the chaotic type and its reheating 
temperature is higher than 10 7 GeV, then the inflationary stochastic gravitational wave 
background would be detected by future space-based interferometric detectors like DECIGO 
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or BBO. Unfortunately we have a small chance to determine the exact value of the reheating 
temperature from the detection. The detection only allows us to set a lower bound on the 
reheating temperature. However, this lower bound 10 7 GeV provides useful and unique 
information of the very early Universe. In this section, as an example, we consider the 
implications of the lower bound of the reheating temperature for the nature of gravitino 
production in the early Universe. 

Many models of supersymmetry breaking, in the context of either supergravity or su- 
perstring theory, predict the presence of scalar fields with Planck-suppressed couplings and 
masses around or heavier than the weak scale. These fields are generically called moduli. 
The coherent oscillation of the modulus soon dominates the Universe, and the late decay of 
the modulus results in very low reheating temperature, upsetting the success of the big-bang 
nucleosynthesis (BBN) [called the moduli problem 39[] unless the modulus is ultraheavy: 
M x > 100 TeV. 

Gravitino is the fermionic superpartner of graviton and has Planck-suppressed interaction. 
The gravitino, once produced, decays with a very long lifetime if it is unstable. The gravitino 
may cause several problems in cosmology. For example, if the gravitino is unstable (if its 
mass is M 3 / 2 ~ lOOGeV — lOOTeV), the decay products of the gravitino would destroy 
the primordial light elements by photodissociation and hadrodissociation, and thus spoil 
the success of BBN (called the gravitino problem 40j). Hence, the yield of the gravitinos 
Y3/2 = n 3/2/s (s is the entropy density) should be constrained: Y3/2 < K5/2 , where the 
value of ^3/™ depends on the gravitino mass. According to the recent analysis 4]J, l42| . 
Kj BN ~ 10" 16 for M 3/2 ~ ITeV and Kg BN ~ 10~ 15 - 10~ 13 for M 3/2 ~ 10 - 100 TeV. On the 
other hand, if the gravitino is stable (for M 3 / 2 < lGeV), the gravitinos can be cold/warm 
dark matter. Hence the abundance of the gravitinos f2 3 /2/i 2 should also be bounded in order 
not to exceed the dark matter abundance fldmh 2 — 0.134 which is precisely determined by 



the WMAP 



Recently, it has been 



ings at the re. 
and moduli) 



leating 



bund that gravitinos are produced not only by the thermal scatter- 



43j but also by the decay of heavy scalar fields (for example, inflaton 



441 ] . Therefore, the "ultraheavy moduli solution" to the moduli problem may 



cause instead a new gravitino problem by the moduli decay. 

In the following, we consider the implication of the possible lower bound of the reheating 
temperature on gravitino cosmology in the context of chaotic/natural inflation. Then the 
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nonthermal production of the gravitinos is only due to the moduli decay 45[ since the 



vacuum expectation value of the inflaton is vanishing for chaotic/natural inflation (with Z 2 
symmetry) 4f|. A similar consideration (without moduli decay) is given in 12]. 



Unstable Gravitino 



Firstly we consider the case of unstable gravitinos. The yield of the gravitinos produced 
by the thermal scatterings at the reheating temperature Trh is estimated by 4l|, |43| 4 

Trr 



Kj* ~ 1.4 x 1(T 12 



10 10 GeV 



(26) 



Thus if the reheating temperature is found to be high (T RH ^ 10 7 GeV), it would immediately 
imply that large gravitino mass M3/2 > 10 TeV is favored 42]. 



The yield of the gravitinos by the X decay is given by 

1 



44 



45 1 



Y x , ~ 

J 3/2 — 



90 



d\ l2 M\ 



192tt V 7r 2 g*(T x ) T x M Pl 



(27) 



where M P1 = mpi/\/87r is the reduced Planck mass and d 3 / 2 — |(X)|/M P1 is related to the 
partial decay rate of the process X — > ip 3 / 2 + ^3/2 as 



3/2 



^3/2 Mj 
288?r M| x ' 



(28) 



Tx is the reheating temperature by X decay with the decay rate Tx — M|/(87rMp 1 ), 5 and 
is given by 



T 



x 



/ QQ \ !/ 4 f M \ 

U^y) ^« = 5 - 8xlo4GeV (ioSilv) 



3/2 



(29) 



where g*(Tx) is the effective relativistic degrees of freedom at T = Tx and we use g*(T > 
ITeV) ~ 220. Then Eq. (J27J) becomes 



x 
3/2 



2.4 x 10 



M 



x 



10 10 GeV 



1/2 



d 



3/2 



(30) 



4 The thermal production of gravitinos at Tx is a factor of px / P<f> smaller than that at Trh and hence may 
be negligible. 

5 We hereby have fixed the order one coefficient. 
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Therefore, from the success of the BBN (Y£ 2 < YE 2 ®), the upper bound on the moduli 
mass is found: 6 

fy BBN \ 2 
Wis) ■ ( 31 ) 



B. Stable Gravitino 



Next, we consider the case of stable gravitinos. The present-day abundance of the grav 



itinos produced by the thermal scatterings at the reheating is given by [43] 



^-^(h^v)^)" 1 ' (32) 

where we have set the gluino mass Mg = ITeV. From Eq. ( )30|) . the abundance of the 
gravitinos produced by the X decay is given by 

h 2 - M ^ ,? ~ 6 8 x l^d 2 , ( Mx ) (^-\ (33) 

where p cr is the present critical density of the Universe and So is the present entropy density. 
The total abundance should satisfy the bound Vt^^h 2 = QjJ^h 2 + Qy 2 h 2 < Qdmh 2 — 0.134. 

Therefore, if Trh > Trh fr° m gravitational wave experiments, then, using Eq. ( )32l) from 
Q™h 2 < fldmh 2 , we find the lower bound on the gravitino mass: 



M 3/2 > 0.21GeV -^f- . (34) 

61 Vl0 7 GeV/ V 0.134 J K ' 

Moreover, from the geometric mean, 

/ / T \ 1/2 / M \ 1/4 

n dm h 2 > n 3/2 h 2 = njf 2 h 2 + n* 2 h 2 > 2^™^ * s.e \w4ev) 

we obtain the upper bound on the moduli mass: 

M X < 6 x 10 'GeV (^) (^)^- (36) 



6 Here we have neglected the effect of dilution by possible entropy productions by moduli decay. The effect 
would decrease the thermal yield by the dilution factor and would weaken the bound on Mx by F 2 

121- 
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V. SUMMARY 



The direct detection of the inflationary gravitational wave background by the next gen- 
eration space-based gravitational wave missions may enable us to explore the early Universe 
more deeply than current observations. In this paper, the detectability of the inflationary 
gravitational wave background in the future experiments, DECIGO and BBO, has been 
estimated using our precise predictions for the spectrum. We have considered several in- 
flation models and have taken into account the effect of the reheating temperature which 
determines the frequency where the signature of reheating arises on the gravitational wave 
spectrum. 

We have found that DECIGO could detect the inflationary gravitational background 
with SNR> 3, for the chaotic inflation model with T RH > 10 7 GeV. The higher sensitivity 
of BBO makes possible a detection with SNR> 5, for the chaotic inflation model with 
Tru > 2 x 10 6 GeV or the natural inflation model with / > mpi and Trh > 10 7 GeV. This 
means, conversely, that non detection of the inflationary gravitational wave background by 
DECIGO would exclude the chaotic inflation model unless the reheating temperature is 
lower than 10 7 GeV. BBO would exclude it unless T RH <2x 10 6 GeV and further exclude 
the natural inflation model with / > mp\ unless T RH < 10 7 GeV. 

We have also discussed the implications of the possible lower bound of the reheating 
temperature on gravitino cosmology in the context of chaotic/natural inflation. Taking into 
account of both thermal and nonthermal production of gravitinos, we find that from the 
lower bound on the reheating temperature we could obtain a lower bound on the gravitino 
mass and an upper bound on the moduli mass. These bounds may provide information 
regarding the gravitino mass and the moduli mass complementary to collider experiments. 
Thus, in future, the direct detection of the inflationary gravitational wave background could 
play a crucial role in probing not only the history of the early Universe but also particle 
physics. 
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